Spatial Structure of Spin Polarons in the t-J Model 
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The deformation of the quantum Neel state induced by a spin polaron is analyzed in a slave fermion 
approach. Our method is based on the selfconsistent Born approximation for Green's and the 
wave function for the quasiparticle. The results of various spin-correlation functions relative to the 
position of the moving hole are discussed and shown to agree with those available from small cluster 
calculations. Antiferromagnetic correlations in the direct neighborhood of the hole are reduced, but 
they remain antiferromagnetic even for J as small as O.li . These correlation functions exhibit dipolar 
distortions in the spin structure, which sensitively depend on the momentum of the quasiparticle. 
Their asymptotic decay with the distance from the hole is governed by power laws, yet the spectral 
weight of the quasiparticles does not vanish. 



PACS numbers; 71.10.Fd, 71.10.-w, 71.27.-ha, 74.90. -l-n 



I. INTRODUCTION 

The problem of spin polarons moving in a quantum an- 
tiferromagnet has found considerable attention, since it 
is important for the description of Mott insulators at low 
doping [|I|. While the major part of investigations for the 
t-J model was concerned e.g. with the polaron dispersion 
and the spectral function using a variety of techniques 
such as exact diagonalization |g-0], selfconsistent Born 
approximation (SCB A) ||-|ll|, string theory ||l|,|l|l and 
other methods ||l4|--|l7| , — our focus here is on the spatial 
structure of the spin polarization and its asymptotic be- 
haviour. The study of the deformation of the spin system 
due to spin polaron formation was mainly performed by 
exact diagonalization techniques plJlq]. However there 
are important questions which can only be studied by an- 
alytical approaches, such as the asymptotic decay of the 
polarization of the medium [|l9|,^ . The latter property 
is closely related to the question whether a quasiparticle 
(QP) description appKes. The first successful measure- 
ment of the single hole dispersion in the Mott insulator 
Sr2Cu02Cl2 by angular resolved photoemission Q has 
revived this interest, and stimulated investigations of the 
t-t'-J model [p3| and more complex Hamiltonians | |23[ |. 

The Green's function for a hole moving in a fixed spin 
background was discussed already in the context of tran- 
sition metal oxides in the late 60th by Bulaevskii et al. 
p3 and by Brinkman and Rice p^ . In those approaches 
the Green's function turned out local and fully inco- 
herent. The first prediction that the low-energy single 
particle excitations in the 2D t — J model l^^l ^i^d its 
anisotropic generalization {Q < a <1), 
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are propagating quasiparticles (QP) with a bandwidth of 
order J was made by Kane, Lee and Read & and was 
confirmed by a number of exact diagonalization studies 
[g|,||. The problem is complicated due to the constraint 
on the fermion operators c] ^ — cj ^(l — rii^-a-) and by the 
fact that quantum fluctuations play a crucial role. This 
model has been widely studied particularly because it is 
believed to contain much of the low-energy physics of the 
high- Tc superconductors p^,p|. 

Nevertheless fundamental issues are still unclear such 
as the spin-dynamics and the form of the Fermi surface 
at moderate doping, i.e. in the regime corresponding to 
the underdoped high-temperature superconductors. But 
even in the case of a single hole there are different views 
e.g. whether the quasiparticle spectral weight is finite or 
vanishes in the thermodynamic limit. In particular An- 
derson has argued that holes introduce a deformation in 
the spin-background which decays as a power law and as 
a consequence the spectral weight should vanish, — lead- 
ing to non- Fermi liquid behaviour p^ . According to this 
argument the non-Fermi liquid behaviour is connected 
with the property of a single hole. Recently Weng et al. 
|28| argued that the quasiparticle weight Zk should van- 
ish as a consequence of string formation associated with 
the Marshall's sign, which is a characteristic property of 
the undoped Heisenberg ground state. These arguments 
are based on the appearance of an orthogonality catas- 
trophe in the matrix element (^k^'"^*|ckCT|0), between the 
exact, i.e. fully relaxed, single hole ground state and the 
state Ckcr|0), where |0) is the ground state of the Heisen- 
berg model without holes. 

The asymptotic decay of the polarization cloud cannot 
be analyzed by numerical methods, such as exact diago- 
nalization (quantum Monte Carlo results for the 2D t-J 
model are still not available) , since such studies are con- 
fined to small clusters and thus can only provide insight 
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into the short-range deformation of the spin-background. 

A particularly powerful tool in the study of the spin 
polaron problem is the slave fermion approach combined 
with a selfconsistent Born approximation for the calcula- 
tion of the polaron Green's function [||J|]. This approach 
was successful in reproducing the diagonalization results 
for the full Green's function obtained by diagonalization 
B. Therefore we shall follow this route here. Further- 
more the method properly accounts for the low-energy 
spin excitations, which are crucial for the long-range dis- 
tortion of the spin-background around the moving hole. 
This method has been also applied to the finite doping 
case 1 29 -3^]. A further important step was the explicit 
construction of the quasiparticle wave function within 
the SCBA by Reiter l^^l- This wave function contains 
implicitely all information about the deformation of the 
spin system, and can be used to calculate this perturba- 
tion in terms of correlation functions. 

Of particular interest is here the study of relative cor- 
relation functions (RCF), i.e. relative to the position of 
the hole, like for example Cr = {no (Sr^ -Sr^)), which 
measures the nearest neighbor correlation function for a 
bond at a distance R = (Ri-fR2)/2 from the hole at 
R = (assuming here that Ri and R2 differ by a lattice 
unit vector u) . Such correlation functions are usually not 
studied because of their complexity. However they pro- 
vide detailed information about deformation of the spin 
system around the moving hole, in contrast to the aver- 
aged correlation function (Srj-Srj), which measures only 
the global change in spin correlations due to the holes. 

The results for the RCF's clearly show that the nearest 
neighbor spin correlations in the neighborhood of the hole 
are reduced, yet they remain antiferromagnetic (even for 
J as small as O.lt). Therefore the frequently invoked 
ferromagnetic polaron picture, where the hole is assumed 
to move in a ferromagnetically aligned neighborhood of 
spins, does not apply to the t-J model. 

The main purpose of this work is to use Reiter 's wave 
function for the calculation of correlation functions and 
to present a quantitative picture of the shape and size 
of the quasiparticle. While a short summary of selected 
results was given earlier [g^ , the present work focuses on 
the description of the technique employed for the calcula- 
tion of the correlation functions. The technique discussed 
here may also be useful in other cases where the non- 
crossing approximation is employed, such as more com- 
plex models including electron-phonon coupling p3| , p4| . 
Results for various correlation functions describing the 
deformation of the spin-background around the hole will 
be presented for the t-J model [a = 1 in Eq. (1)] as well as 
for the simpler t-J^ (a — 0) model |35| which has no spin- 
dynamics and has a simple classical Neel ground state. 
For the t-J model the relative correlation functions are 
found to be strongly dependent on the momentum of the 
quasiparticle and in good agreement with known results 
from exact diagonalization. 



Furthermore a detailed investigation of the asymp- 
totic decay of various correlation functions is given. For 
example the perturbation of the nearest-neighbor spin- 
correlation function Cr is found to decay as 1/R'^ with 
the distance from the hole. Since the asymptotic be- 
haviour of these correlation functions is closely connected 
with the question whether Z^ is finite or not, it is im- 
portant to calculate the deformation of the spin-system 
within the different existing approaches. In the present 
framework it is found that all perturbations introduced 
by the hole in the quantum antiferromagnet decay at 
large distance as power-law with dipolar or more com- 
plex angular dependence depending on the momentum 
of the quasiparticle. Nevertheless this does not lead to 
vanishing quasiparticle spectral weight, consistent with 
earlier numerical results based on the study of the po- 
laron Green's function within the SCBA ||lO(| . 

The plan of the paper is as follows: In Sections II and 
III we briefly summarise the selfconsistent Born treat- 
ment for the Green's and the wave function of the quasi- 
particle, and provide the framework to calculate expec- 
tation values with respect to Reiter's wave function. Sec- 
tion IV deals with the quasiparticle spectral weight, the 
magnon distribution function and provides a discussion 
of the convergency of the approach. The more complex 
RCF's are studied in Section V for two generic cases, 
the t-J and the t-J^ model, i.e. one with and the other 
without spin dynamics. This section also contains a dis- 
cussion of the asymptotic behaviour of the different cor- 
relation functions. The paper concludes with a summary 
in Section VI. 



II. SLAVE FERMION APPROACH 

In a first step of the reformulation of the problem, holes 
are described as spinless (slave) fermion operators, i.e. 
on the A-sublattice a spinless fermion creation operator 
is defined as hf — Ci-\ while the corresponding operator 
Cii — hf Sf is expressed as a composite operator, and 
similarly for the _B-sublattice [Q. The kinetic energy 
then consists of terms of form —thjhtSj, that is, each 
hop of the fermion is connected with a spin-flip. The 
spin dynamics is described within linear spin wave theory 
(LSW) which provides a satisfactory approximation for 
the 2D spin-1/2 Heisenberg antiferromagnet. 

We follow here Refs. ||^- [lO| , |36|| and express spin op- 
erators via the Holstein-Primakoff transformation, and 
simplify the notation by performing a 180° rotation of 
the spins on the S-sublatticc, 
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Here the origin Rq = belongs to ^-sublattice (spin 
up) and Q = {n / a,T: /a). The lattice constant is a = 1. 
The spin interaction term is further diagonalized after 
linearizing spin operators and performing the Bogoliubov 
transformation 



u,q 
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where foj, — N^^^^J^i^^^'^'^'^H ^'iid N is the number 
of lattice sites. Here we use the usual Bogoliubov co- 
efficients 



"-qj 



and the spin wave dispersion is cjq — 
2Jy/l — (a7q)2 with 7q — {cosq^ + cosqy)/2. After 
fermion operators are decoupled into slave fermions and 
bosons, 
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the fermion-magnon Hamiltonian emerges |^ 
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A constant term irrelevant for the present discussion has 
been dropped here. One recognises that the kinetic en- 
ergy appears now as a fermion-magnon coupling with a 
coupling function given by Afkq = 4i(uq7k-q + Wq7k)- 
This Hamiltonian is similar to the small polaron model 
except that a kinetic energy term for the spinless fermions 
is absent. In the case of the cuprate superconductors, 
where t > J , the model is in the intermediate or strong 
coupling regime and a selfconsistent calculation tech- 
nique must therefore be chosen. 

In the following we will use the hole Green's function 

where the QP band energy ek and the pole strength Zk 
are related to the fermion self-energy SkC"^) as ek = 
Sk(ek) and Z^^ = 1 — dT,i^{uj)/duj\f^ respectively. 

We calculate I]k(u;) within the self-consistent Born ap- 
proximation (SCBA) 
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Such an approximation amounts to the summation of 
non-crossing diagrams to all orders. The validity of this 
approach is well established. The QP dispersion and 
spectral weight calculated within the SCBA [Q agrees 
very well with the exact diagonalization results for small 
cluster s p ]. The spectral weight in the limit iV — > oo is 
finite [|10{. In the extreme J/t '3> 1 limit |M, however, 
this method leads to Zk — > 1, i.e. an overestimation in 
comparison to 0.82 obtained for t = in Ref. p^. The 
success of the SCBA has roots in the vanishing of low or- 
der vertex corrections as pointed out by several authors 
for systems where the hole is coupled to an AFM spin 
background [p6|P,||. 



III. QUASIPARTICLE WAVE FUNCTION 

Given the Green's function in selfconsistent Born ap- 
proximation it would be interesting to have the wave 
function of the quasiparticle which corresponds to the 
pole in Eq. (^) at energy ek = Sk(ek)- The knowledge 
of this wave function will allow us to calculate in princi- 
ple all equal-time correlation functions which define the 
perturbation of the AFM-background around the hole. 
We follow here closely Reiter's |Q original approach and 
prove in addition that the quasiparticle weight derived 
from the wave function is consistent with the well known 
expression obtained from the Green's function. 

The quasiparticle wave function is defined as the eigen- 
state of H 
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which gives rise to the quasiparticle peak in the spectral 
representation for the Green's function 
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Here |0) represents the vacuum state with respect to 
fermion and magnon operators and |^km) is an eigen- 
state of Hamiltonian Eq. (@) with eigenenergy ekm. The 
spectral weight of the quasiparticle state |^k) 

^k -I (*k \h+\0) p (10) 

can be quite small, however it should not scale to zero in 
the thermodynamic limit, whereas the matrix elements 
contributing to the incoherent part are of 0(1/A^) or 
smaller. 

Given the Hamiltonian Eq. (ph we expect the quasi- 
particle wave function | ^k) to have the form 

I vPk) = a°(k)/.+ |0) + ^Eai(k,qi)/i^q,a+ |0) 
"^ qi 

— Y^ a2(k,qi,q2)/ik-qi-q2aq2aqj0) 
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+ .. 



(11) 



where the coefficients a"(k, qi, ..., q„) are to be deter- 
mined. 

/^From the Schrodinger equation we obtain the follow- 
ing system of equations for the expansion coefficients: 
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To obtain these equations which correspond to the non- 
crossing approximation for the Green's function one 
has to adopt the following contraction rule: When one 
magnon is annihilated in the n-magnon component of the 
wave function, Eq. (11), only the contribution is consid- 
ered where the last magnon in the sequence, i.e. a't , is 
annihilated. This is reminescent of the retraceable path 
approximation in momentum space. The general equa- 
tion for n > reads; 

- Ljqi ... -tjq,Ja"(k, ...,q„) - a""^(k, ..., q„)Mk„_i,q„ 

- ^ E a"+i(k,...,q„+i)Mk„,q„+, =0, (14) 

q™+i 

where k„ = k — qi — . . . — q„ . 

As first shown by Reiter [621 this sequence of equations 
(O)-(O) has the general solution 
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a"+^(k, ...,q„+i) = a"(k, ...,q„)gk„,q„+i , 
where we introduced the abbreviation 

5k„.q„ + i = Mk„,q„ + iGk„+i(t^ - CJq, - ... - OJq^+J ■ (16) 

Substituting Eq. ( |l5| ) into the last term on the l.h.s. 
of Eq. (|lj), we recognise that this term is identical to 
the expression Eq. (H) for the selfenergy Sk('-^) times a". 
This yields for Eq. (nj) 
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,qn-i)Mk„_i,q„=0. (17) 



Since the prefactor of a" is the inverse of the Green's 



function G^^ (w — Wq^ - 



this equation is identical 



to Eq. (|15| ) with n replaced hy n — 1. It only remains to 
be shown that also Eq. (^2|) is solved. Equation (|lj) 
becomes 



a°(k)(a'-Ek(w)) =0, 



(18) 



which has a nontrivial solution a"(k) 7^ at the QP- 
energy lo = e^. The knowledge of the Green's function 
Eq. (^ is sufficient to calculate from Eq. ( p^ ) iteratively 
the coefficients a"(k, qi, ...q„). 

The coefficient a"(k) which determines the QP-weight 
Z]i = [a"(k)]^ follows from the normalisation of the wave 
function (vl/k | ^k) = ^^=0 I «"(k, -.qn) P= 1, 

(*k|*k)-[a°(k)](l + l^ 
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When one calculates the derivative dT,]^{uj) / duj from 
Eq. (Q) and compares the result with Eq. (n9h it is easy 
to see that \M 
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As I^Pk) is normalized to 1, [a"(k)]'^ is indeed identical 
to the QP-spectral weight as calculated directly from G. 
This latter step is important, since it accomplishes the 
prove of the internal consistency of G and vp, i.e. where 
both are calculated within selfconsistent Born approxi- 
mation. It should be emphasised that the above deriva- 
tion does not rely on the assumption that the coupling 
term in the Hamiltonian is small. 

Because of the presence of AFM long-range order the 
quasiparticles move on one sublattice, while visiting the 
other sublattice only virtually. In view of the 'degener- 
acy' Ek+Q = Ck and Gk+Q = Gk we define new Bloch- 
operators 
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which create holes on the t (i) sublattice for r = 1(— 1), 
respectively. The momenta k are now restricted to the 
reduced (AFM) Brillouin zone. The corresponding wave 
functions including magnon operators up to order n are 
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Here gk.q — AfkqGk-q(ek — f^^q) as defined in Eq. (|lq). 
We note that the Green's function G in gk,q is always 
evaluated below the lowest pole and therefore real. For 
example in Gk-q(ek — Wq) the energy ek-q > Ck — i^q, 
hence 5k,q is real. This is actually a subtle consequence 
of the selfconsistent evaluation of the Green's function 
which leads to a smaller energy variation of the QP en- 
ergy compared to the spin wave dispersion |IQ] . We 
stress that this holds also true in the strong coupling case 
where also the ek variation is of order J, i.e. comparable 
with spin wave energies. The choice of sublattice wave 
functions Eq. (02) is convenient since they are eigenstates 



of 5fot with eigenvalues ±1/2. 

The diagrammatic structure of the wavefunction l^^ki ) 
is shown in Fig. 1(a). The translation rules are straight- 
forward: (1) Open ends on the right correspond to oper- 
ators h] and aj, ...aj, , (2) thin lines are as- 

k— qi — ... — q„ qi q^i ' ^ -' 

sociatcd with V^, (3) a vertex (dot) connected with a 
double line corresponds to gk.q, and (4) there is a mo- 
mentum sum for each magnon line. It is obvious that 



the wavefunction does not correspond to a strict order n 
expansion with respect to the fermion-magnon coupUng, 
since the Green's function, Fig. 1(b), is already evaluated 
selfconsistently with respect to this interaction. 

In the next Section we will investigate the relative 
importance of the different terms in the wave function 
Eq. (pS) and address the question under which condi- 
tions this series can be truncated. 



at J^ /t ~ 0.3. For smaller J^/i the number of magnon 
terms needed to fulfil the sum rule A/k = 1 increases 
rapidly. In Fig. 4 the distribution of magnons A^ is 
displayed for the strong coupling case, J^/f <C 1 . In this 
regime A^ has a maximum at a finite value n, which 
increases with the coupling constant t/J^. 

The average number of magnons forming the spin po- 
laron may be defined as 



IV. MAGNON DISTRIBUTION FUNCTION 
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The first question one may ask is: "How many 
magnons are involved in the formation of the polaron?" 
As the coupling between hole and spin-excitations is the 
kinetic energy of the t-J model, small values of J/t cor- 
respond to strong coupling (small spin stiffness), where 
many magnons are excited by the hole motion. In order 
to estimate the number n of magnon terms needed in the 
wave function we have calculated the norm A/k 
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The distribution function A^^ defines the probability 
for the n-magnon contribution in the wave function. A 
similar study on a small cluster was presented in Ref. Q . 
In Fig. 2 the norm is presented diagrammatically con- 
sistent with Eq. ([l9|). Each term A^ corresponds to a 
single non-crossing diagram with n-magnons. Vertices 
denoted with dots correspond to the fermion-magnon 
coupling matrix elements Mkq and the double line to 
the square of Green's function Eq. (p) calculated within 
SCBA. The analytical expression for A^' is independent 
of T and given by 
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for m > 0, while A\^ = Z\^. From Eq. ( |19[ ) we know that 
A/k — > 1 in the limit n — *■ oo. This normalisation condi- 
tion will serve as a check of our numerical procedure. 

It is instructive to study the distribution function A-^ 
first in the case of the the t-J^' model (a = 0). In this 
limit of model Eq. (H^) the analysis becomes simple be- 
cause there is no intrinsic spin-dynamics. The SCBA 
equations for the self-energy are independent of k and 
reduce to one equation Sk(w) — 4t^[cj — 2J^ — I]k(^ ^ 
2J^)]^^ H]. Equation (BS) can then be expressed in a 
recurrence form 
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The norm A/k is shown in Fig. 3 as a function of the 
number of magnon terms n for various J^ /t. A crossover 
between the weak and the strong coupling regime occurs 



In the Ising limit (n) is identical to the average number of 
spin deviations (local magnons) (^- Sf S~) = (^^ a\ai). 
It is evident that the latter expression is proportional to 
the average string length lav |19|1 of overturned spins in 
the Neel state created by the hole motion. As the string 
potential is an approximately linear function of the string 
length this implies (n) ex lav « {t/ J'^Y^'^. This estimate 
is reasonable for long strings, i.e. small J^/t. In Fig. 5 we 
present (n) as a function of J^ /t calculated with up to 40 
magnon terms in the wave function. For large J^ /t^ 1 
only the leading term to = 1 in Eq. ( P7| ) is relevant, 
therefore the asymptotic result is {n) = (i/J^)^. For 
J^/t <^ 1 we find excellent agreement with the result 
(n) = 1.4(t/J^)i/3 obtained by Mattis and Chen ||l|- 

From these results for the t-J'^ model it is clear for 
J^/t > 0.4 the wave function can be truncated at n = 3 
or even at n = 2. We note that the same holds true for 
the i- J model ^. 

In Figures 6(a) and 6(b) the numerical results for the 
norm A/k of the t-J model are shown for k = and 
k = (-1, -1), respectively, both calculated with up to n = 
3 magnons kept in the wavefunction. For J = 0.4 3- 
niagnon contributions are necessary to fulfil the norm. 
The quasiparticle spectral weight Zk, which corresponds 
to the n = term displayed in Figs. 6(a,b), is always 
finite except in the limit of vanishing spin stiffness J — 0. 
Thus our wavefunction docs not lead to an orthogonality 
catastrophe. This result will be further complemented in 
the next Section by a detailed study of the asymptotic 
decay of the spin-polaron correlation functions. 

The question whether the QP spectral weight Zk for 
the t-J model is finite or not is still not completely set- 
tled. Numerical results obtained on small clusters are in a 
good agreement with the results obtained from the SCBA 
Equations d^ (lO). In the SCBA formahsm Eq. (|), Zk is 
finite |P,pq,nO| because the holc-magnon coupling matrix 
element for g — + is not singular and therefore the hole is 
weakly coupled to low energy spin waves. In Ref. pq | it 
was argued that Zk should vanish nevertheless because of 
string like phases associated with the hole motion (due to 
hidden Marshall signs). We stress that the Mar hall sign 
convention is implicitely included in our present formu- 
lation. In fact the vacuum state |0) (originating via uni- 
tary transformation) is equivalent to the quantum Neel 



groundstate of the T = Heisenberg model, and thus by 
construction obeys the Marshall sign rule in the original 
basis, i.e. before the 180° rotation of the B-sublattice. 
After the transformation, Eq. (g), there are no additional 
phases in the transformed Hamiltonian due to the Mar- 
shall sign. 

Not considered in the present treatment, is the effect 
of the 4 broken bonds meeting at the site of the hole. 
This leads to an additional relaxation of the spin cor- 
relations and hence to a reduction of the quasiparticle 
weight. This effect is expected to be strongest in the 
limit t = Q. The exact result for the spectral weight in 
this case is Z = 0.82 and was derived by Mal'shukov 
and Mahan jSg] (as compared to 1 in the present treat- 
ment). The energy change due to the broken bonds must 
also be included in the Born approximation if one wants 
to compare the quasiparticle energies with those from ex- 
act diagonalization, as discussed by Martinez and Horsch 



V. SPIN-POLARON CORRELATION 
FUNCTIONS 

The spatial structure of the spin-polaron can be de- 
scribed with various correlation functions measuring the 
perturbation of the spin system relative to the position 
of the moving hole. As we shall see, these correlation 
functions are strikingly different in the t-J and the t-J^ 
model, — a consequence of the absence of spin dynamics 
in the latter model. In the t-J model perturbations cre- 
ated by the hole are carried away by spin waves thereby 
generating a power law perturbation pattern with an in- 
teresting angular dependence, whereas in the absence of 
spin dynamics the perturbations are characterized by an 
isotropic gaussian decay. 

Such relative correlation functions can be evaluated 
using the quasiparticle wave function. One of the sim- 
plest correlation functions is the distribution of magnons 
around the hole 
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(no(a^aR)). 
(28) 



Here n.i — h\hi is density operator for holes at site i 
with position R.;. A'r, also corresponds to the distribu- 
tion of spin deviations, {no{S^S^)) . Therefore it pro- 
vides a suitable measure of the polaron size. This cor- 
relation function is also proportional to the distributions 

Correlation functions such as A^r are evaluated using 
similar diagrams as in the calculation of the norm A/k. 
One has to evaluate the expectation values 
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Here the summation ^^ corresponds to all lattice sites 
and the density operator for the hole 



kik2 



(30) 



has to be expressed in terms of operators h\^r , Eq- (121 
The operator Or, is decomposed into magnon variables 



^^^ ^Yl [/qiq2(R-)"qi"q2+5qiq2(R-)aqi"q2- 

qiq2 



-H.c. 



The diagrammatic structure of the contributions for a 
general correlation function of this type is presented in 
Fig. 7. The first class of diagrams is symmetric and de- 
rives from the vertex function /q^qj (R). These diagrams, 
denoted by (-Bn), arise as diagonal contributions from the 
n-magnon component of the wave function. 

The construction rule for these diagrams is the follow- 
ing: If the vertex / (circle) as well as the connected two 
magnon lines (together with their vertices and associated 
double lines) are removed from the diagram, one must 
arrive at a diagram contained in the expression for the 
norm (Fig. 2). Otherwise the diagram is not consistent 
with the selfconsistent Born approximation and should 
be dropped. 

The second class of diagrams {Cnm) is asymmetric and 
corresponds to the vertex function g^-^^^^^R) which con- 
nects n magnon contributions with tti = n ± 2 magnon 
terms in the wave function. Again only such diagrams 
must be taken into account which are consistent with 
the construction rule formulated before. 

The vertex functions / and g are expressed in terms 
of Bogoliubov coefficients and thus strongly momentum 
dependent. For the case of the correlation function A'r 
we have 

/q,q,(R) = i(tiq,«q, +Z;q,Z;qJe'('l^-1=)-^ 

5qiq,(R) - i(Uq,t;q, +«q,UqJe'('l^+'l^)-^. (31) 

In order to illustrate a typical calculation of matrix el- 
ements needed in the correlation functions, we present 
here the second order contributions B2 in Fig. 7, 

B2 = N~^ 2Z /qiq2(R-)(3k,q3gk-q3,qi5k-q3,q2ffk,q3 + 



qiq2q3 

+ 3k,qi5k-qi,q33k— q2,q35k,q2 ji 



(32) 



where the first and the second term correspond to non- 
crossing and crossing term, respectively. 



A. Ising limit (a = 0) 

In general correlation functions and the corresponding 
matrix elements have to be evaluated numerically, which 



is easy for not too large n up to ~ 5. The t-J^ model is an 
exception, since the Bogoliubov factors simplify to Uq = 
1 and -Wq = 0. Thus Re /q^q2(R) — cos(qi — q2) ■ R 
and (7qiq2(R) = 0, respectively. As the Green's func- 
tion is dispersionless it is possible to express the matrix 
elements analytically and perform the summation of di- 
agrams {Bn) to any order. Furthermore diagrams {Cnm) 
are zero. It is instructive to express the wave function 
in real space. Each n-magnon term can then be visu- 
alized as a string of n steps with starting point at the 
origin. From such a study one can get insight into the 
noncrossing structure of the wave function and correla- 
tion functions. 

The SCBA is similar to the retraceable path approxi- 
mation, yet with the important difference that in SCBA 
the hole can also hop backwards on its path. At the level 
of the Green's function the differences were discussed in 
Rcf. |l^. The result for the magnon distribution func- 
tion, Eq. (Ea) can therefore be expressed as 



m— 1 



,0) 



mPn 



(33) 



where Pm — X]?=m^k^ '"^^ ^^ interpreted as a prob- 
ability to have at least m local magnons excited. The 
coefficients pm (R) represent the probability that a string 
of m excited local magnons ends at a given lattice posi- 
tion R. This distribution can be determined by counting 
all possible paths of m steps, where in each step all z 
neighbors can be reached. 



Pm(R) = 4- 



m 



(34) 



Here m± — {m — \ R^ \ ± | Ry | ) /2 must be a non- negative 
integer, otherwise Pm(R) = 0. This result is free of 
boundary conditions. 

The correlation function N^ can be used to deter- 
mine the spatial size of the polaron in the Ising limit. 
We define the size of the polaron quantitatively by the 
radius Rp (element of Bravais lattice), which encloses 
a given fraction p of the total number of spin devia- 
tions, p — {n)^^J2R<R ^R.- I'^ Fig. 8 the polaron 
radius Rp vs. J/t is shown for three different values 
of p = 0.75, 0.9, and 0.99. In the physically inter- 
esting regime, J^/t ^ 0.3, the polaron is contained 



\l/2 



oc 



within the radius R < 2. The scaling Rp oc (n 
[t/J'^Y^^ expected for the polaron [E9pfl] is well es- 
tablished. We have also calculated the average radius, 
{R) — (?^)^^ X]r l-f^l-^R-j ^^"i ttie root-mean square ra- 
dius, Rrms = ((n)-iER|R|'^R)'/'- In Fig. 8 Rrms 
and {R) are presented with solid and dashed lines, re- 
spectively. The RMS radius can be well fitted with 
R = 1.06(i/ J)"i" for J^/t < 1. 

In the Ising limit the total spin is not conserved. How- 
ever, the z-component of spin is a conserved quantity. A 



state with one static hole {t = 0) at the site io has by def- 



inition the z-component of total spin S^^ 
— t/2 (t = ±1), i.e. the spin of one site of the sublat- 
tice not corresponding to io- If the hole becomes mobile 
(t ^ 0), some spins around the hole deviate from the Neel 
order. The region where the spin order is disturbed cor- 
responds to the spin polaron defined above. The corre- 
lation function describing the spatial distribution of spin 
around the hole is thus 



5R 



,Q R y^ ^ 



re 



iQR; 



(^i«"R.+R"R.>+R.) ' (35) 



where we have expressed spin operators in terms of 
magnons according to Eq. (0). The conservation of spin 
corresponds to the sum rule X]i?,7to ^^ ~ "''"Z^- "Fhc local 
spin operator is within the LSW approximation related 
to the number of bosons, Sf = -^ — ajai. However, iSr 
is due to the factor exp(zQ • R^) non-trivially related to 
TVr and has to be calculated independently. After carry- 
ing out the steps similar as in the evaluation of A^r one 
obtains 



5, 



R 



JQR 



9^0 



p„i(R)P„ 



(36) 



andP,„ = E;=™(-l)^"4'^ 

The spin correlation function S'r for several J^/t val- 
ues is given in Fig. 9. We have performed the calculation 
for n up to 40, which was more than sufficient to obtain 
converged values. The results can be qualitatively un- 
derstood visualising the correlation function 5r in the 
moving coordinate frame of the hole. For large J^/i the 
hole moves slowly through the Neel ordered background 
and on the average spends more time on sublattice r. 
The alternating contribution to S'r corresponds to the 
AFM ordered background, and is given by the first term 
on the r.h.s of Eq. (|3^), which is apart from the AF- 
alternation independent of R. It represents the differ- 
ence in the probability that the hole sits on the t and 
I sublattice, respectively. This background contribution 
tends to zero for J^/i ^ 1, where the hole rapidly hops 
from one sublattice to the other. The second term in S'r 
carries all spatial dependence, i.e. defines the region of 
spin disturbance, and becomes dominant at J^ /t <^ 1. 

We would like to stress here that the disappearance of 
the staggered Neel structure for small J^ /t in this corre- 
lation function is simply a consequence of the fact that 
the hole visits the two sublattices with equal probability, 
and it does not mean that the antiferromagnetic order 
is no longer present as one could naively conclude from 
similar results of a finite cluster diagonalizations. We 
note that, our results resemble surprisingly well the re- 
sults for S'r obtained in exact diagonalization studies for 
small clusters f0,Ef . 



In Fig. 10 we show with open squares the dependence 
of Sr at R = (1,1) with J^/i = 0.4 on the num- 
ber of magnons n taken into account in the calculation. 
The results for other J^/t values are in agreement with 
the results for A/k, where we found that above (below) 
J^ /t ^ 0.3 a relatively small (large) number of magnons 
are excited and therefore needed in the evaluation of the 
correlation functions. 

The conservation of the total spin z-component can be 
tested by summing up X]i?=£o "^R-- '^'^^ total spin S^^^ is 
presented in Fig. 11 as a function of J^/i with diamonds 
and the full line is a guide to the eye. S^^^ consists of two 



parts. The first corresponds to the first term in Eq. (36) 



\Po, and is shown with the dashed curve. Pq represents 
the difference in the probability of the hole sitting on 
sublattice t or |. The second term in Eq. ( |3q ) is not 
presented separately. The interchange of importance of 
the two contributions is in agreement with the discussion 
above. The small violation of the S^^i conservation law is 
a consequence of the Holstein-Primakoff representation of 
spin operators. We have also calculated S^^^ as a function 
of n. For J^/t > 0.3 only three magnon terms included 
in the wave function give sufficient accuracy in agreement 
with calculation of the norm A/k- 



B. Heisenberg limit (a = 1) 

The important new features of the t-J model are (i) 
the spin-dynamics described by antiferromagnetic spin 
waves, which have a linear dispersion around q = (0, 0) 
and (tt, 7r), respectively, (ii) The ground state of the 
model in 2D is a quantum Neel state, i.e. more com- 
plex than the simple classical Neel ground state of the 
t-J^ model. An immediate consequence of (i) is that a 
spin-deviation which is created by a single move of the 
hole will propagate away from the hole in form of a spin- 
wave until it is reabsorbed at a later instance. The long 
wavelength spin excitations determine the distortion of 
the quantum antifcrromagnet at large distances from the 
hole. 

A further aim of our study of the RCF's is to show 
that the spin correlations remain antiferromagnetic in 
the vicinity of the hole. The antiferromagnetic corre- 
lations are weakened yet not ferromagnetic. The ferro- 
magnetic polaron picture, i.e. a carrier accompanied by 
a ferromagnetically aligned spin-cluster, does not apply 
here. Ferromagnetic polarons are a quite popular sce- 
nario usually inferred by a generalization of Nagaoka's 
theorem 03] , which applies to the J = model, to finite 
exchange interaction J. 

To gain more insight into the complex angular depen- 
dence of the relative correlation functions we present in 
addition to the numerical results (which include up to 
n=3 magnons) an analytical study of the RCF's based 
on the wavefunction in the one-magnon approximation. 



This wave function is sufficient for a quantitative discus- 
sion in the large J case; yet it also predicts the large 
distance behaviour for smaller J values. 

The main k-dependence in the wave function stems 
from the hole-magnon coupling matrix element Mkq. In 
the q —t Q limit the k- and q- dependence of gik.q oc 
g^^/^(7k — 2^/^Vk • q/?) determines the asymptotic sym- 
metry of the correlation functions. From this structure 
of .gk.q it is clear that at fc = the spatial symmetry 
is s-wave, whereas at the minimum of the QP band at 
k = (-S, -f) the symmetry is determined by the dipolar 
term, where Vk =Vk7k (S^Q. 

If one is only concerned about the behaviour of the 
wave function at large distance R from the position of 
the hole R.; the one magnon contribution simplifies and 
one can perform the corresponding Fourier transform of 
.9k, q- The Bloch representation of the wave function in 
the limits t/ J —* and R ^ oo in then approximated in 
leading order. 
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Here the Fourier transforms (jio — —2^/2^\^t / {J R) and 
(/>i = — 2(vk •R)i/(Ji?^) have different spatial symme- 
tries. The (/)i-term is dipolar and vanishes at k = (0, 0) 
and (7r,7r). At (-I,-!) 4>i has its maximum, while the 
monopole contribution 0o vanishes instead. We note that 
|^jj.| ) has similarity to the wave function describing the 
motion of a '^He atom in superfluid ''He [^,Q. In the 
following this wave function will serve us as a starting 
point for the derivation of the asymptotic properties of 
various correlation functions. 

The wave function, Eq. (p3|), is properly normalized 
also for the Heisenberg limit and the norm is given by 
Eqs. (H) and (H). The evaluation of A^'^ can be done 
numerically. In Fig. 12 A^ is plotted for k — (fc, k) and 

n = 0, 1,2,3 at J/t = 0.4. For n = 0, A^ = Zk and 
the momentum dependence is well known pC| . The next 
term, n = 1, corresponds to the emission of one magnon. 
The momentum dependence is very weak, which can be 
qualitatively understood from the t/ J —> limit. For 
q < qc -^ 1 the one-magnon contribution A\^ follows as: 



A[ 



U^^ 



9l. 



qgd(7(x(72 + |vk|')gc. (38) 



Here we have put Zk ^ 1 for the weak coupling limit. The 
obtained result is constant for k along the (1,1) line. This 
behavior is found in the full numerical calculations even 
in the strong coupling regime J/t = 0.4 in Fig. 12. Other 
distribution functions A^^ in Fig. 12 have a more subtle 



momentum dependence which cannot be reproduced with 
this simple asymptotic expansion. The sum X]"=o ^k 
is close to 1, as it is clear also from Figs. 6(a,b). The 
results in Fig. 12 show that the higher order magnon 
terms are less important for quasiparticle momenta close 
to the band minimum at k — (^, ^). For the full J/t 
dependence of the norm at k = (^, ^) see Ref. [|o|. In 
order to obtain converged results in the Heisenberg limit, 
we have performed numerical calculations using unit cells 
with A^ = 16 X 16 up to A^ = 32 X 32. In summations 
over the Brillouin zone the points q = and q = Q were 
excluded. The numerical method of solving the SCBA 
equations for G'k(w) was identical to Ref. [ p6[ . 

The average number of magnons (n) in Fig. 13 is pre- 
sented for J/t = 0.4 and momentum k — (^, ^), i.e. 
close to the QP band minimum. It is interesting that {n) 
calculated for the t-J model almost coincides with the 
result obtained for the Ising case (Fig. 5). 

The additional spin-deviations created by the hole mo- 
tion are given by the expression A^r = (no(a^aR)) — 
TVafm- Here we have subtracted the large contribution 
Nafm = 0.197 due to the quantum fluctuations in the 
ground state of the 2D Heisenberg antiferromagnet in 
the absence of the hole. The shape of the polaron is 
elongated in the direction of the QP momentum which 
reflects a quasi one-dimensional motion of the polaron, as 
was pointed out in Ref. [gO[. This is consistent with the 
asymmetry of the QP energy band in the "hole pocket" 
region centered around k = (f , f ), where the effective 
next-nearest neighbor hopping for the (1,1) direction is 
^ 5x that in the (1,-1) direction. This asymmetry is 
most pronounced at the bottom of the QP band and grad- 
ually vanishes away from the QP energy minimum and 
disappears at A: = and k = (tt, 0). In the limit i? — > oo 
the perturbative result is to lowest order in t/ J given by 
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,vk R, 



(39) 



This result strictly holds only asymptotically, but nev- 
ertheless it reflects all symmetries found at short dis- 
tances in the numerical treatment. The momentum de- 
pendence is qualitatively correct as well, while the J/t 
dependence is correct only in the range of validity of per- 
turbation, t/J -^ 0. The correlation function A^r decays 
as a power-law, A'r cx R^^. Although the number of ex- 
cited magnons (n) is small, it turns out that the change 
in the total number of spin deviations ^^ N-r diverges 
logarithmically. The definition of the polaron size used 
for the Ising limit of the model thus cannot be used here. 
Since the magnetic excitations cjq vanish linearly with q, 
also other correlation functions show power law decay, 
yet with different exponents EQ]. 

In Fig. 14 we display the distribution of z component 
of spin S'r as a function of J/t and for k = (^i ^)- 
This correlation function depends strongly on the direc- 
tion and size of the momentum of the quasiparticle (see 



Ref. iQ ) ■ The asymmetry of the polaron is reflected in 
different values for S'r at positions R labelled with 2 and 
2' or 4 and 4', respectively. This result is quite different 
from the isotropic perturbation in the Ising limit (e.g.. 
Fig. 9). We stress that the same asymmetry was found 
in numerical studies of an 18 site t-J cluster with one hole 
||l8|. The ground state is at k = (|7r, fvr) for J/t = 0.4 
[The point k = (f , f ) is absent in that system]. Due to 
the high symmetry of the 4x4 cluster such subtle asym- 
metries of the polaron cannot be studied there. The J/t 
dependence of S'r is in excellent agreement with the re- 
sults of Refs. [0,0. As finite size effects in such small 
clusters are expected to be quite large, the good agree- 
ment of Sr with the exact results is surprising. The rea- 
son for the disappearance of the AFM structure in this 
correlation function for small J/t is as in the Ising limit a 
consequence of fast hole motion which leads to an average 
over the two sublattices. It does not imply that the anti- 
ferromagnetic order of the spin background is destroyed. 
The correlation function is small in this limit because the 
polaron is large and many sites contribute to the sum rule 

qz 1 

'-'tot — 2' 

To test our analytical and numerical procedure we have 
calculated Sr for different number of magnons terms n 
(Fig. 10). The convergence rate is similar as in the case of 
Ising limit, i.e. three magnon terms give a sufficient ac- 
curacy. To display the anisotropy Sr is shown in Fig. 10 
for R = (1,±1) and k = (f,f) with J/t = 0.4 as a 
function of number of magnon lines in the wave function, 
n — 1,2, 3. The corresponding contributions (diagrams) 
are labelled with symbols (-Bn) and (Cnm) as defined in 
Fig. 7. The asymmetry of the polaron, which is fully 
consistent with the numerical results of Ref. [Q , can be 
attributed to the diagram C02, corresponding to a two 
magnon process. 

The effect of the hole on the AF correlations and the 
energy of the spin system is measured by the nearest- 
neighbor spin correlation function Cr — {uq (Sr^ -Srj)) 
defined on bonds between two neighboring sites (1 — 2), 
R = (Ri + R2)/2 ||. In Fig. 15(a) Cr is shown as a 
function of J/t and for k = (^,^)- The correlation 
function remains negative and in agreement with the nu- 
merical result obtained on a 16 sites cluster [g[. Hence 
AFM-correlations persist in the vicinity of the hole con- 
trary to what one would expect from the ferromagnetic 
polaron picture. Cr is asymmetric as can be seen, e.g., 
from the bonds R — (1,±|), labeled with 1 and 1'. The 
momentum dependence of Cr can be explained with the 
perturbative result which follows from asymptotic wave 
function Eq. (|^) M 
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(40) 



This correlation function decays as R~'^ at large distances 
and could be used as a definition of the size of the po- 
laron. Our results suggest that the size of the polaron 



measured with this correlation function is, at moderate 
J /t = 0.4, of the order of a few lattice sites. 

Another interesting aspect of the deformation of the 
spin background is contained in the bond-spin currents 
JR = ("o(Sri X Sr2)^u), where u is a unit vector u — 
R2 — Ri [|l8|,Q . This quantity follows from the equation 
of motion for the spin density 



>R 



: it'^{&ss' Ck^^cr+u.s' - H.c.) - 2i,/^ Sr x Sr+u, 



where a are Pauli spin matrices. Here the first term is 
the spin current induced by the hopping of the hole and 
the second term (~ j/j) describes the backflow in the spin 
system. Due to the broken symmetry total spin is not a 
good quantum number, therefore we consider only the z 
component of the current. In Fig. 15(b) jr is presented 
as a function of J/t and for k = (^, ^). Jr is an odd 
function with respect to the wave vector (at fc = 0). Be- 
cause of symmetry it vanishes also at k = (7i',0). Since 
the ground state has AFM long-range order, the points k 
and k -I- (tt, tt) are equivalent, and therefore, Jr vanishes 
also at k = {^, ^). Comparison of jr with exact diag- 
onalization results is delicate. As reported in Ref. |2y] 
we find good agreement with results from Ref. ||l^. For 
the complete momentum dependence of Jr in compari- 
son with exact diagonalization see Ref. Q|. Agreement 
is excellent in the anisotropic limit, a < 1. In the Heisen- 
berg limit a reliable comparison is very difficult because 
of the strong k dependence of Jr which makes it very 
sensitive to the boundary conditions of small clusters. In 
Fig. 15(b) we present jr for various bonds R defined in 
Fig. 15(a). The asymptotic pattern of bond spin currents 
is dipolar [pO|, 
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■7k [Vk - 



(vk-R)R . 
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(41) 



The spin backfiow current Jr decays as R~^ and vanishes 
in the ground state for k = {^,^). Again the general 
momentum dependence is correct as in the case of the 
other correlation functions considered. 



VI. CONCLUSIONS 

We have outlined a method that allows to calculate 
the real-space structure of a spin polaron in a quan- 
tum antiferromagnet. The approach is based on the 
spin-polaron formulation for the t-J model, where holes 
are described as spinlcss fcrmions while the spin exci- 
tations are treated within linear spin-wave theory. The 
single-particle Green's function in (selfconsistcnt) Born 
approximation, which has been shown earlier to provide 
an excellent description of the numerical data for the t- 
J model, is sufficient to calculate the many-body wave 
function describing the polaron (Reiter's wave function). 



We have shown here how this wave function can be used 
to calculate in the frame of the selfconsistcnt Born ap- 
proximation quite complex correlation functions. 

Our calculation of a number of correlation functions, 
which measure the deformation of the spin system due 
to a moving hole (spin polaron) , provides a very detailed 
check of this approach against exact diagonalization. In 
particular we have shown that the spectral weight of the 
spin-polaron quasiparticle calculated from the wave func- 
tion is consistent with the result derived from the Green's 
function. It is demonstrated how the number of spin- 
excitations involved in the polaron formation increases 
with decreasing spin stiffness when J — > 0. We have de- 
termined the probability distribution of the number of 
magnons excited in the ground state and found that for 
J/t ~ 0.4 (a typical value for copper-oxide superconduc- 
tors) the average number of magnons is about one. In 
the Ising limit the average number of magnons scales as 
(n) oc {t/ J^Y^'^ as J^/t — > in agreement with the string 
picture for the moving hole and with the Ref. |{4l|| . 

We have put particular emphasis on the asymptotic de- 
cay of the perturbations introduced by the spin polaron 
formation. Since the spin wave energies Wq in a quantum 
antiferromagnet vanish linearly with g, perturbations in 
the spin system decay with a power law. For example 
the change of the local spin deviations TVr cx i?~^, while 
the perturbation of the nearest-neighbor spin correlations 
decays as R~'^ with the distance from the hole. In the 
t-J model all correlation functions have a quite complex 
structure in real space which depends on the momentum 
of the polaron, whereas in the i- J^ model all perturba- 
tions are isotropic. 

We note that despite of the power law decay of the 
polaron correlation functions in the t-J model the quasi- 
particle spectral weight does not vanish. Whether this is 
correct or an artifact of the selfconsistcnt Born approx- 
imation remains to be shown by a more rigorous treat- 
ment. 

Finally we want to stress that the approach discussed 
here may also be applied to other interesting problems 
such as strongly correlated electrons coupled to Holstein 
or other phonons. 
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FIG. 1. Diagrammatic representation of the wave function 
l*k" ) '^'-"" ''' ~T- (a-) The first three terms contain no-magnon, 
one-magnon and two-magnons excitations, respectively, (b) 
The double line represents the single particle Green's function 
in noncrossing approximation. 



FIG. 2. Diagrammatic representation of the norm TVk. The 
first term - containing no magnon line - is identical to the 
QP pole strength Zk. 



FIG. 3. The norm TVk in the Ising limit as a function of 
the number, n, of magnons included in the wave function for 
various J^/t. 



(u) 



FIG. 4. The distribution of the number of magnons A^^ 
as a function of n for various J^/t. 



FIG. 5. Average number of magnons {n) in the Ising limit 
as a function of J^/t. The inset shows (n) on logarithmic 
scale compared with asymptotic results. 
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FIG. 6. The norm Ak for the t-J model as a function of 
J/t including magnons up to order n = 3 in the quasiparticle 
wavefunction. (a) Momentum at the top of the QP band at 
k = 0, and (b) at k = (^, ^) , i.e. close to the QP-band 
minimum. The numerical calculation was performed for a 
momentum space grid corresponding to a 32 x 32 system. 
The solid and the dashed lines are guides to the eye only. 



FIG. 7. Diagrammatic representation of correlation func- 
tions. Each class of {B„) diagrams contains n noncrossing 
diagrams. (C„m) diagrams appear always in pairs with the 
H.c. counterparts. For a detailed description see text. 



FIG. 14. The 0-component of the spin correlation function 
5r vs. J/t for k = (^, ^). Note the asymmetry between 
the directions R||(l, 1) and R||(l, —1) in the t-J model. The 
numerical calculation was performed on a grid corresponding 
to 16 X 16 system. 



FIG. 15. (a) Nearest-neighbor spin-correlation function Cr 
and (b) the ^-component of the bond spin currents Jr as func- 
tion of J/t for the quasiparticle momentum k = (%-,%-). The 
inset in (a) provides a definition of the n.n.-correlations con- 
sidered. In both (a) and (b) note the asymmetry between the 
directions R||(l, 1) and R||(l, — 1). The asymptotic behavior 
is given by Eqs. (40) and (41), respectively. 



FIG. 8. Radius of the polaron in the Ising limit vs. J^/t 
for various definitions: Rp for p — 0.75 (full circles), p — 0.90 
(open squares), and full squares for p = 0.99. Here p defines 
the fraction of spin deviations within the radius Rp. The 
full line represents root mean square radius, Rrms, while the 
dashed line gives the average radius (R). 



FIG. 9. Distribution of the 2;-component of spin 
Sn = (woSr) around the moving hole for three different val- 
ues for J^/f. 



FIG. 10. The dependence of Sr on the number of magnon 
terms n in the wave function. In the Ising limit ^(i^i) 
(dash-dotted line) is essentially converged for n > 3 given 
J" /t = 0.4. In the Heisenberg case the contributions from 
different diagrams to S'r are shown for R = (1,1) (solid) 
and R = (1,-1) (dashed line), respectively, for J/t — 0.4 
and k = (|-, |-). For the classification of diagrams (Bn) and 
(C„m) see Fig. 7. 



FIG. 11. The total ^-component of spin Stot vs. J" /t (di- 
amonds). The dashed line represents iPo in Eq. (36), where 
Po is the difference of the probabilities for the hole to occupy 
sublattice | or J,, respectively. 



FIG. 12. Magnon distribution function ^j^" for the t-J 
model {J/t = 0.4) as function of k = (fc, k) for various magnon 
numbers n: n — (diamonds) n — 1, (full circles), n — 2 
(open squares), and n = 3 (open circles). The numerical cal- 
culation of all matrix elements was performed on a grid cor- 
responding to a 16 X 16 system. Lines connecting the symbols 
are guides to the eye only. 



FIG. 13. The average number of magnons (n) involved in 
the spin polaron formation in the t-J model as function of 
J/t. The polaron momentum is k = (^, ^). 
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